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ABSTRACT 


Optical Matrix Vector Multipliers (OMVM) are of interest to real time adaptive 
array beam forming. Feasibility of high speed compact signal processing for 
applications such as sidelobe cancelling has yet to be demonstrated. The 
importance of the OMVM stems from the following: 

(1) Because of its optical parallelism it can operate at extremely high 
speeds, surpassing by several orders of magnitude competing electronic 
techniques . 

(2) The matrix-vector multiplication operation is basic to a variety of 
important signal processing applications, including adaptive array 
processing and other general control operations. 

The very high speed of the OMVM strongly suggests its application to iterative 
processing, where in a very short time the system can cycle through perhaps 100 
iterations, and assiiming necessary convergence, yield a solution. This method 
has been suggested as a means of inverting the covariance matrix associated 
with noise samples from adaptive array antenna elements. A review of the jammer 
problem for radars indicates that solution is required in the 1 ysec to 10 ]isec 
range. The spatial light modulator must be updated within this time frame. 

We will describe here the optical architecture of an incoherent matrix vector 
multiplier where the matrix is complex and separable into outer product vectors . 
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INTRODUCTION 


In the late 1970’ s, researchers at Stanford University began investigation 
of a new class of incoherent electro-optical processors for high speed matrix- 
vector multipliers (ref. 1, 2). Work has continued, both at Stanford (ref. 3) 
and Carnegie Mellon Institute, to the point were significant processing 
capabilities have been demonstrated. It is not unreasonable to expect that 
such systems could perform 10^ matrix-vector multiplies per second with 100x100 
element matrices. 

One limitation of these processors as developed thus far is the relatively 
fixed nature of the mask representing the matrix. in general photographic 
masks have been employed, and processing where the matrix must change with some 
rapidity has been impossible. 

In our search for a matrix mask that can be updated rapidly, we have considered 
several candidate spatial light modulators (ref. 4) (SLMs) , including the 
Itek PROM, and the Hughes LCLV. Neither the PROM nor the LCLV can be cycled 
at rates approaching 1 MHz, and these candidates are thus unacceptable. A mosaic 
of pockels cells could perhaps be made to operate in the desired fashion. 

However, there are the usual problems of switching high voltages in microseconds, 
and severe electronic crosstalk appears inevitable. "We 

have settled on a scheme whereby the mask is implemented acousto-optically . 

Acousto-optic modulators possess several desirable features when used in signal 
processing systems. These include reliability, availability, high response 
speeds , and wide bandwidths . However , they are not easily adaptable to use as 2-D 
spatial light modulators. The basic technique and system architecture for a 
pair of 1-D SLMs is discussed in the following section. 


ACOUSTO-OPTIC IMPLEMENTATION OF VECTOR MULTIPLICATION BY OUTER PRODUCT MATRICES 


The operation to be implemented is: 


Y =: H X (1) 

where X is a complex input column vector. 


H is a complex matrix separable into outer products: 

H = A 


( 2 ) 


The multiplication of separable matrices would be possible if the OMVM were a 
coherent processor. When the necessary complex to non-negative real encoding 
of H is implemented for an incoherent processor, the separable form of the system 
matrix is destroyed. However by using time multiplexing we can still avoid the 
requirement for 2-D masks. Consider one possible method. We can define 
expressions in which vectors and matrices are defined in the following form (ref 5,6). 
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where X^ , X^, X^; Y^, Y^, Y^; Hq, H 2 are all non-negative 

real and where 0 = 2 7T/3. The components X^, Y^, H^, i = 0,1,2 have the 
dimensions of X, Y, and H, respectively. 


The basic matrix vector multiplication operation can be written: 
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Consider the case where H is an outer product matrix (2) . If we express 

- ^ ^ 

(2) in terms of a three component decomposition, we have H = (A^ + A^e + (5) 
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Multiplying through and collecting terms In powers of exp (i9) we find that 
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To represent H by a three-component decomposition, we need nine outer 
products. 
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These nine outer products can be produced using a pair of crossed acousto - 
optic cells with multiple transducers. Figure 1 illustrates the case where 
A is a two-dimensional complex matrix. The desired outer products are 

formed at the intersection of the acoustic waves, in the following format: 
Let M = 


0 0 0 1 


M = 


^0 


^ «0 '^2 "l 


\ ^2 


\ ^2 


(7) 


2 2 


At this stage we encounter a difficulty in system architecture: although_the 

various outer products of (3), together with input vector components X^, 
X^ f sufficient to allow calculation of Y, these outer products are not 

arranged so as to allow direct implementation of equation (4) . The three- 
component form of matrix H 
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is itself not representable as an outer product. As a 

consequence, it appears as though system architecture will of necessity 
be more complicated than would be the case were H directly realized by, 
for example, a more general 2-D SLM. 

Let us consider the problem further to gain a better understanding of its 
nature. We simplify notation by denoting the elements of equation (7) by 
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From equations (3) , (6 ) , (7) and (9) we have 
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Inspection shows that component can be expressed by 


Y = 
n 


I 

i, j fk 


M. ,X, 
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where i/j,k satisfy the relationship 

(i+j+k) mod 3 = n 
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Given the matrix M and input X = (X^ possible to generate 

all the terms of (10) by a single vector -matrix multiplication operation. 

It is, however, possible to generate Y with three vector-matrix multiplication 
operations , as follows . 

First, input the vector (X X calculate its product with M. 

We write the resulting product as 
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Y, = M X^ 
la — 2 


'2a 


This is followed by the operations 
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If the outputs of these three operations are properly combined, the result is 
the desired three-component representation for Y. Specifically, 
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study of these equations shows that as the components r and X^ are 

stepped cyclically in the input, the components of Y^, Y^, and Y^ are 

themselves stepped cyclically. One method for performing the simmiation of 
(12) and (13) is to use charge -coup led integrating detectors in the output 
plane of the optical system. The integrated charges can then be stepped 
synchronously with the input components. Another method is to use high 
speed analog switches and short delay lines for storage. 


PLANS FOR IMPLEMENTATION 


The basic optical layout of the iterative processor is shown schematically 
in figure 2. Light emitted from the sources is collimated into N parallel 
beams by means of a set lens Ll. Each lens is a segment of a spherical lens. 

Thus, each row of the first SLM is illuminated by a collimated beam from one 
of the light sources. The output beams from the AO Bragg cell, A, are focused 
by lens L4, spatially filtered by SI, and recollimated by L5. The lens 
set L4-L5 images the aperture of Bragg cell A at the second Bragg cell, B. 

The modulated beam from the second Bragg cell is spatially filtered and brought 
to focus on the linear detector array. The spatial fillers are necessary 
only to remove the zero and second order beams. The output beams are modulated 
in intensity at NxN points within the optical aperture. This representation 
is equivalent to having both masks located at the same plane. The outputs 
from the light emitters are modulated^by vector^s whose elements are X . . On 

multiplication, the output terms are Y. = V** A. B. X.. ^ 

1 2^ D 1 D 

j 

When a non-negative real approach is implemented, the above scheme must be 
modified. The speed of the processor is slowed down by approximately a factor 
of three, since for each iteration of three successive multiplications must be 
carried out. After each multiplication the analog data must be stored while the 
input vector is shifted. The three successive outputs are added as in equation 
(18) . The output vector then becomes the input to the light emitters as in other 
iterative optical processors . 

The time multiplexed method described can be implemented using CCD shift registers. 
Unfortunately, state of the art CCD shift registers operating at the necessary 
speed lack the dynamic range and noise properties necessary for o\ir applications. 

We have thus chosen to implement the previously mentioned procedure utilizing fast 
four channel CMOS analog multiplexer switches arranged in a switching matrix con- 
figuration. Referring to figure 3, as the MUX is clocked, its output is shifted 
along the linear laser digde arrays. Similarly, at the linear pj^otodiode array, 
the outputs ^2a^ from the first input cycle (X^ X^ X^) 

are stored^in coaxial cable delay lines. ^ The input X vector is then switched to 
(X X X^) and the output (Y^^^ Y^^ ^Ob^ switched to the proper location 

and stored in another delay line. This procedure is repeated a third time, 
with the appropriate combinations taken to yield the desired outputs. 
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Some of the problems particular to the implementation involve non-linearities 
in the Bragg diffraction process, electrical and acoustic crosstalk, and limited 
dynamic range due to light scattering. 


CONCLUSIONS 


An analysis has been performed of the applicability of iterative optical 
processors in environments requiring a short transient response time. An 
approach has been proposed utilizing realizable acousto-optic single 
dimension light modulators where the matrix is the outer product of two vectors. 
The problems associated with a non-negative real representation have been 
addressed. A solution has been proposed albeit at some cost of processor 
speed. A proof of concept experiment is underway which will better determine 
the limitations and advantages of this method. 
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N COLLIMATED BEAMS BEAM REDUCER MASK FORMED BY N CHANNEL 

laUMINATING ROWS BRAGG CELL MODULATORS AND SPATIAL FILTERS 
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Figure 2.- Optical system for matrix vector 
processor with separable masks. 
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